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Abstract 

Static topologically-nontrivial configurations in sigma-models, for spatial dimension 
D >2, are unstable. The question addressed here is whether such sigma-model soli- 
tons can be stabilized by steady rotation in internal space; that is, rotation in a global 
S0(2) symmetry. This is the mechanism which stabilizes Q-balls (non-topological 
solitons). The conclusion is that the Q-mechanism can stabilize topological solitons 
in D = 2 spatial dimensions, but not for D = 3. 
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1 Introduction 



Topological solitons (stable, localized, topologically non-trivial solutions in a field 
theory) have long been of great interest, both for their mathematical properties and 
for their applications in many areas of physics. In any system admitting such solitons, 
the non-trivial topology is not sufficient to ensure that the solitons are stable. One 
obvious reason for this is that topology is, by definition, a non-metric structure, and 
so it cannot determine the size of the solitons; for that, one needs to balance the 
forces acting on the soliton in such a way that it has a preferred size. Recall, for 
example, the 0(n) sigma model in D-l-l dimensions (with trivial boundary condition 
at spatial infinity), which is the subject of this article. This system admits topological 
configurations (textures) whenever the homotopy group tt£){S^^^) is non-trivial; in 
particular, for {D,n) equal to (2,3), (3,3) or (3,4). If L) = 3, then solitons tend to 
shrink — in the pure sigma model, there are no static solutions. In the D = 2 case, 
there are static solutions, for example the Belavin-Polyakov solitons [Ij in the 0(3) 
system; but these are unstable [21 El- 

A soliton can always be prevented (or rather discouraged) from spreading out by 
the addition, if necessary, of a potential (a term involving only the field, and not its 
gradient). In order to stabilize the soliton size, we also need to introduce something 
which prevents it from shrinking. There are several possibilities for such an anti- 
shrinking mechanism: for example, a Skyrme term involving four (or more) powers 
of the field gradient; or a gauge field suitably coupled to the sigma-field; or periodic 
time- dependence (rotation in an internal space). This third possibility also underlies 
non-top ological solitons (Q-balls). In this paper, we investigate to what extent the 
Q-ball mechanism is effective at stabilizing topological sigma-model solitons. We 
shall see that stationary topological Q-solitons exist in D = 2 spatial dimensions, 
but not for D = 3. This result is analogous to that for the Landau-Lifshitz equation 
(Heisenberg model of ferromagnetism) , as one might have surmised since the static 
Landau-Lifshitz system is identical to the static sigma model. 



The 0(3) sigma model involves a scalar field taking values on 5^; this field can be 
represented as a unit 3- vector cp = (0i, (j)2,4'3) with (j) ■ (j) = 1. The Lagrangian is 



where V is some potential function (which, for simplicity, we take to depend only 
on (p^). The space-time coordinates are = {t, x^), with j = 1, . . . , D. The system 



2 Q-Solitons in the 0(3) Sigma Model 





2 



has a global S0(3) symmetry which is broken, by the potential term, to S0(2). This 
SO (2) acts only on (pi and 02, namely by changing the phase oi (p := 4>i + i(/)2- The 
corresponding conserved quantity is 

Q = j lm(00)(i^x. (2) 

Minimizing the energy of a configuration subject to Q being fixed implies 0j, in 
particular, that (p has the form 

cP{t,x^)=i;ix^)e'''\ <^3 = V.3(x^'), (3) 

with I^P + (V's)^ = 1- Without loss of generality, we shall assume that > 0. Note 
that Q = 0/, where / = / IV'P d^x. The energy of a configuration of the form Q is 
E = Ed + Ek + Ep, where 




The boundary condition is V'3 ^ 1 as r — > oo; so we need V'(l) = 0. 

A stationary Q-lump is a critical point of the energy functional subject 
to Q having some fixed value. Such a Q-lump is (classically) stable if this critical 
point is a local minimum of E. The usual (Derrick) scaling argument shows that any 
stationary Q-lump must satisfy 

(2 - D)Ea - DEp + DE^ = 0. (4) 

Let the positive constant m be defined by V'{\) = — m^; in other words, ^^(V's) ~ 
m^(l — V's) ~ |m^|'i/'P for ~ 1. Then, near spatial infinity, the Euler-Lagrange 
equations corresponding to E imply that 

- [m^ - 0^)V' = 0. 

So in order to satisfy the boundary condition ip ^ Q as r ^ oo, we need < m. The 
solitons are exponentially localized if < m, but less-localized solitons with = m 
may also exist. 

The parameter m is part of the specification of the system, and the parameter Q 
is set by the initial data. Each of these two parameters has dimensions; the combi- 
nation Qm^~^ is dimensionless, whereas the combination (Q/my^^ has dimension 
of length, and determines the size of the soliton. Configurations in this system are 
classified topologically by their topological charge (an integer); ii D = 2, then N 
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is the winding number, while if D = 3, then N is the Hopf number. Let E{N, Q) 
denote the energy of a configuration (or rather, of data </), (fit) with topological charge 
N and Noether charge Q. 

Let Vo be the normalized potential function Vq := ^V/m?. Note that Vo(^3)/|^P — 
1 as ~^ 1- It is clear from @ that if Vq > everywhere, with Vq > IV'P 
somewhere, then there can be no solution. So the constant K defined hy K = 
min [Vo(V'3)/|''/'P] should satisfy K <1. It then follows that 

Ep = 'n? j Vb > rT?KI = {Km^/0^)Ek > {K'n?/0^)Ep, (5) 

where the final inequality comes from (jlj). As a consequence, we have 

Krn^ < 0^ < m^. (6) 

In the = 3 case, the first inequality is strict: Km? < 0^ < m?. 

In two spatial dimensions, it is possible to have K = \^ which corresponds to the 
choice V{ip'^) = \m?'{\ — ijj^). So here = m. This system '5^., and generalizations 
in which the target space is some other Kahler manifold, arise naturally by dimen- 
sional reduction from 'pure' sigma models in one dimension higher IB]. The energy 
satisfies a Bogomolny bound E{N, Q) > 4ttN + mQ, which (for N > 2) can be sat- 
urated: for each value of > 2 and Q, there is an explicit family of stationary 
multi-soliton solutions such that E{N, Q) = 4ttN + mQ. There is no force between 
the individual solitons: in particular, the total energy has the additive property 
E{Ni, Qi) +E{N2, Q2) = E{Ni + N2, Qi + Q2)- One may use moduli-space methods 
(as was done in [3 IH1 El for other sigma-model systems) to investigate the scattering 
of moving solitons [21 ; this involves finite-dimensional mechanics on M . The dynam- 
ics turns out to be rather exotic (as is also the case jlUl 1111 112j for non-topological 
Q-balls). The solitons are only polynomially localized, and the non-existence of an 
N = 1 soliton is related to this; an A = 1 configuration tends to shrink in size, and 
there is no stationary N = 1 solution. 

On the other hand, if AT < 1, then 1-solitons can exist. Different choices of 
V (having K < 1) seem to lead to similar behaviour, but this has yet to be fully 
investigated; in what follows, we take Vq = i(l — V'3), so K = ^. Let us consider, first, 
the thin- wall limit iT^ . where Qm^^^ ^1. In this limit, the (bulk) contributions 
Ep and E^ to the energy are very much greater than the (surface) contribution E^. 
So the energy is approximately E f {1 — Tp^) + ^Q'^ / / (I — V'i)- Without loss 

of generality, we may assume that tp^ = 1 outside of some compact set. So space is 
partitioned into three regions: one (with infinite volume) where ^3 = 1, the second 
(with volume A) where |V'3| 7^ 1, and the third (with volume B) where ip3 = — 1. 
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Note that E depends on A (and on the value that tp^ takes on A), but not on B, and 
that 

— = (02 -mVi)V'3- 
dtps 

So for fixed A, the function E has a minimum for ips = {ips = =b0/m are local 
maxima). Hence we should set ip^ = on A, and the energy becomes 

E = \Q^IA + \Tr?A. 

Thus for a given value of Q, the energy has the minimum value i?min = mQ/^/2 when 
A = Q^/l/m. Note that = m/\/2, at the lower end of its allowed range ®. 




Figure 1: The profile function /(r) for the 1-soliton solution on the plane, with m = 20 
and Q = IGOtt. 



To make further progress, we need to include the effect of surface tension, in 
other words include the term E^. Let us consider, first, the planar case D = 2. For 
simplicity, we assume rotational symmetry about a point in the plane: the field is 
taken to have the form ip = sm{f) exp{iN9) and tp^ = cos(/), where / = /(r). The 
boundary conditions are /(O) = n and /(oo) = 0, and N is the topological charge. 
The energy functional E = E^ + Ej. + Ep was minimized numerically, for various 
values of m, Q and N . The term E). was used in the form E'^ = Q'^ /I, so that one 
can minimize while keeping Q fixed; the quantity lv does not enter explicitly, but can 
be derived (via the formula = Q/I) once the minimum has been found. In each case 
that was investigated, a smooth minimum was reached. For = 1 and Qm = 32007r 
(close to the thin-wall limit), the profile function /(r) is plotted in Fig^ We see that 
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there is a region around r = where /(r) drops rapidly from vr to 7r/2 (the term 
prevents this region from shrinking further in size); and then a region (corresponding 
to A in the argument above) where / = tt/2 4^ ijj^ = 0. Outside of this region, the 
field takes on its asymptotic value f = <^ ip3 = 1. 

Fig 121 displays results for = 1, m = 1 and a range of values of Q. We see that 
E is very close to being linear in Q (recall that in the Bogomolny case, it is exactly 
linear): to a very good approximation, we have E = 4tt + 3Q/4. As for 0, we know 
from © that has to be in the range 1 / \/2 < < 1 , and we see from the figure that 
this is so; furthermore, — > l/\/2 as Q — > cxd (the thin- wall limit |E])) while 0^1 
as Q ^ (the thick- wall limit ^14 J. 
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Figure 2: The energy E and angular frequency of the 1-soliton, as functions of Q. 



In the Bogomolny case jS] mentioned previously, the energy E{N, Q) of a station- 
ary soliton has the feature that, for a given Q, the quantity i?Ar := (47rA^)^^£'(A^, NQ) 
is independent of N: in fact, i?Ar = 1 + mQ/A-K. This corresponds to the fact that 
in this Bogomolny-type system, there is no force between stationary solitons. For 
the potential V = \m?{\ — ^|), however, there are such forces. This can be seen 
by examining E^ for fixed Q = Att and for various values of A^. Tabled shows the 
results for 1 < < 7. The energy density of the 1-soliton is peaked at the point 
r = 0, whereas that for the rotationally-symmetric A^-soliton is peaked on a ring. 
Note that E^ is a decreasing function of N ^ which suggests that the this A^-soliton is 
stable against breakup into solitons of lower topological charge. But this remains to 
be checked; in particular, one should investigate the vibrational modes about these 
rotationally-symmetric solutions . 

Finally, let us turn to the case oi D = 2> spatial dimensions. Configurations with 
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1 2 3 4 5 6 7 



Em 1.888 1.840 1.826 1.820 1.817 1.816 1.815 
0/m 0.859 0.826 0.817 0.814 0.882 0.811 0.810 

Table 1: Energy E]^ = E{N, NQ)/{A'kN) and frequency 0/m, for m = 1 and Q = Air 

nonzero Hopf number N look like closed loops, which may be linked or knotted. 
Static knot-solitons are have been studied in the Faddeev-Skyrme system, where a 
Skyrme term is added to the Lagrangian: this extra term stabilizes the solitons, 
which would otherwise shrink |15| 1161 \17\ I18j . The question here is whether there 
exist stationary Hopf solitons which are stabilized by internal rotation rather than by 
a Skyrme term. The answer to this question appears to be negative; what happens 
is as follows. Consider an = 1 configuration. Typically, -03 = 1 at spatial infinity 
and on a curve which extends to infinity; let us visualize this curve as the z-axis. 
Secondly, ^p3 = —1 on a closed loop L around the z-axis. Finally, "03 = on a torus 
around the z-axis, with the loop L in its interior. So the region A in the previous 
thin-wall analysis is a thickened torus (with the i3-region in its interior) resembling 
a closed string. Numerical experiments indicate that, roughly speaking, the Q-effect 
supports the thickness of the string, but not its length; the string has a tension 
which causes its length to shrink. So the configuration collapses, and there is no 
stationary Hopf soliton in this system. In principle it remains a possibility that for 
some potential function V, some value of Q, and some nonzero value of the Hopf 
number A^, there might exist a stationary solution; but this seems rather unlikely. 

Another way of viewing the situation is as follows. The Q-mechanism provides a 
lower bound on the quantity / (since to is bounded above, and Q = luj is fixed); this 
in turn means that the volume of the soliton is bounded below. But surface tension 
then acts to make the soliton spherical. So we are led to the following conjecture: any 
stationary Q-ball (whether topological or not) with uj < m, in D spatial dimensions, 
has 0(D) symmetry. (In the Bogomolny case, where cu = m, rotational symmetry is 
not essential (Hj.) The instability of Hopf Q-solitons is an immediate consequence of 
this conjecture, since 0(3) symmetry implies that the Hopf number is zero. 

3 The 0(4) Sigma Model in 3+1 Dimensions 

In this section, we investigate the analogous problem for the 0(4) sigma model in 
three space dimensions. The details of the system are similar to those of the previous 
section. The field takes values on S^, and is represented as a unit 4- vector (j) = 
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{(po, 4'i,4'2,4'3)- The Lagrangian is (0), as before, but with the potential being ahowed 
to depend on (pQ and ^3: V = V{cj}Q,(j)3). In general, this breaks the global 0(4) 
symmetry to 0(2), the subgroup which rotates (pi and (f)2, leaving 0o and (p^ fixed. 
The expressions for E^i, E]. and Ep are the same as before, with D = 2>, except that 
in E^ there is an extra term involving (djipo)'^. 

We take the boundary condition to be "00 ~^ 1 as r — > cx) in The mass m is 
defined by V^('0Oi V'a) ~ m'^il — V'o) ^ ^m'^{\tp\'^ + ^'3) for V'o ^ 1; the function Vq and 
the constant K are defined as before. The virial relation from Q with D = 3 holds 
as before, as does the inequality Km? < 0^ < m?. 

In ^H]) this system was studied, with the potential Vq = 2(1 — ^po). Since in that 
case we have K = 1, no soliton solution can exist. The authors of |19j reach this 
conclusion for the topologically trivial case = 0; they report numerical evidence 
for a non-trivial solution with = 1, but this cannot be correct. 

In order to allow the possibility of non-trivial solutions, we need a potential Vq 
which has K < 1; for what follows, we shall take Vq = ^{1 — tPq) as in the previous 
section. Then there are solutions, but it appears that they all have trivial topology 
(A^ = 0). One way to see what happens is to consider the thin- wall limit, where 
m^Q » 1. So the energy is approximately E / (1 — tp^) + / \'>p\'^', and 

the corresponding variational equations are 

CVo-mVo = 0, eV'-0'V' = O, ^^^3 = 0, (7) 

where ^ := mPipQ + 0^\ip>\^ (the ^-term arises from enforcing the constraint (p- (f) = 1). 
These equations ((TJ have a number of solutions, namely: 

• Vo = 0, ^ = 0, Vs = ±1; 
. = 0, V's = 0, Vo = ±1; 
. Vo = V'a = 0, IVP = 1; 

• V3 = 0, ^0 = ±0/m, l^/^p = 1 — 0^ /m?. 

So to construct a mimimum-energy configuration, we must partition space into 
regions (separated by infinitesimally-thin walls), on each of which one of these rela- 
tions holds. It is clear that regions on which V's 7^ contribute only to Ep, and that 
we can reduce the total energy by instead setting V's = 0, V'o = 1 on these regions. 
In other words, ip^ 'collapses' to zero, and is replaced by V'o- 

This is exactly what one sees in numerical simulations. For example, we may 
start with the 0(3)-symmetric 'hedgehog' ansatz 

V'o + iV'j'^j = exp [i/(r) x^Oj/r] ; (8) 
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here aj denotes the Pauh matrices, and the profile function /(r) satisfies the usual 
boundary conditions /(O) = vr, /(oo) = 0. The winding number is = 1. Note 
that tp = tpi + iip2 vanishes on the x^-axis, and that V'o(O) = — 1- If we now relax 
the configuration by flowing down the energy gradient, then ip^ approaches zero 
everywhere except at the single point r = 0; in other words, there is no continuous 
minimum in this topological class. By contrast, there is a smooth minimum which 
has ip3 = 0: this is topologically trivial, and is essentially a standard (nontopological) 
Q-bah. 



4 Concluding Remarks 



We begin with a few remarks on the similarities with stationary topological soliton 
solutions of the Landau-Lifshitz equation 

^ = -<^x— . 9 
dt dc/) 

Here is a unit 3-vector representing the local orientation of magnetization, and the 
energy E is given by 



E 



i(a,0)-(a,</.) + [/(</.3) d^x. (10) 



A typical choice for the function U is U = A{1 — where ^4 is a constant; this 
corresponds to an easy-axis anisotropy. The boundary condition is (p^ ^ 1 as r ^ co. 
The total magnetization 

M = [ (l-(/.3)(i^X (11) 



is a conserved quantity. It is clear from scaling that the only static solutions are 
the Belavin-Polyakov solitons ^ in spatial dimension D = 2, with U = 0. But by 
allowing time dependence, more solutions are possible. In particular, we may allow 
periodic time dependence, and look for stationary solutions such that 

01 + i02 = Hx'W"', h = Mx'). (12) 
With this ansatz, the Landau-Lifshitz equation Q is equivalent to 

-^(E-luM) = 0. (13) 
ocp 

We may think of the solutions as critical points of E, subject to the constraint 
that M has a given value 20 . The simplest example occurs if f/ = |z^(l — ^3), 
for then the functional appearing in ()13() consists of only the gradient term, and 
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the Belavin-Polyakov solitons are solutions (in D = 2); these correspond, of course, 
to the Bogomohiy-type Q-lump sohtons The analysis of stationary topological 
Landau-Lifshitz solitons leads to rather similar results as for Q-solitons (although 
the dynamics of moving solitons is quite different). In D = 2, there are topological 
solutions (called magnetic bubbles — see [U for a review); a single soliton is pinned 
in space, and cannot move i22, . In D = 3, on the other hand, there are no stationary 
Hopf solitons [22]; however, such solitons can be stabilized by allowing them to move 
at constant velocity ^511^ . 

Returning to sigma- model dynamics, we have seen that the Q-mechanism stabi- 
lizes topological solitons in D = 2 spatial dimensions, but not in D = 3. Stabilizing 
vortex rings (Hopf textures) in D = 3 is particularly difficult, since there are two 
length-scales (the length of the loop and its width), each of which has to be fixed. 
The Q-effect can stabilize the latter, but not the former. One does get stable loops 
in systems with a Skyrme term , and also in systems with a magnetic field 

sufficiently strongly coupled to the scalar field (minimal coupling is not enough) |25j . 
But in the basic versions of each of these systems, there is only one length-scale; and 
so the length of the loop is of the same order as (and only slightly greater than) its 
thickness. It remains an open question as to whether there is a system admitting a 
stable Hopf soliton in which the two length-scales are significantly different. 

A sigma-model soliton in D = 2 can be thought of as a (straight) sigma-model 
string in three spatial dimensions. So, for example, the Q-stabilized solitons discussed 
in this paper may find application as cosmic strings. Given an appropriate potential 
V, long strings with internal rotational energy will be stable, although closed loops 
will eventually shrink and decay. In this connection, it is worth recalling that, on a 
cosmological scale, both the width and the length of sigma-model strings are stabi- 
lized by cosmological expansion; but stabilizing the length requires a greater rate of 
expansion than stabilizing the width |26j . 
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